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1. Introduction 

The theoretical description of weakly bound exotic nuclei close to the drip-line is one 
of the most exciting challenges today. What makes this subject both particularly 
interesting and difficult, is the proximity of the particle continuum implying strong 
modification of the effective nucleon-nucleon interaction and causing unusual spatial 
properties of the nucleon density distribution (halo structures, large diffusivity) . Many 
of those nuclei are involved in the chain of thermonuclear reactions and, in the absence 
of data at relevant energies, the models of stars rely to certain extent on calculated 
astrophysical factors (see Bahcall (1989)) . 

In weakly bound exotic systems, the number of excited bound states or narrow 
resonances is small and, moreover, they couple strongly to the particle continuum. 
Hence, these systems should be described in the quantum open system formalism which 
does not artificially separate the subspaces of (quasi-) bound (the Q-subspace) and 
scattering (the P-subspace) states. For well bound nuclei close to the /3-stability line, 
microscopic description of states in the first subspace is given by nuclear SM with model- 
space dependent effective two-body interactions, whereas the latter subspace is treated 
in terms of coupled channels equations. In this work, we question the validity of this 
basic paradigm of nuclear physics, and propose its modification for weakly stable exotic 
nuclei by taking into account coupling between Q and P subspaces in terms of residual 
nucleon-nucleon interaction. This coupling modifies the scattering solutions as well as 
the spectroscopic quantities for interior bound states. 

As said before, we are interested in describing low lying bound and quasi-bound 
states in exotic nuclei. For that reason, we can restrict description of particle continuum 
to the subset of one nucleon decay channels. Still, in few rare cases of two-nucleon halo 
nuclei, this limitation may turn out to be restrictive. In any case, further improvement 
of our model to more complicated channels like, e.g., a - channels, can be done as well 
(Balashov et al (1964)). 

2. The shell model embedded in the continuum (SMEC) 

The influence of scattering continuum on the mean-field properties has been addressed 
before (see Dobaczewski et al (1996) for review). No such analysis has been done so 
far for the realistic SM. A possible starting point could be the Continuum Shell Model 
(CSM) approach (Fano 1961, Mahaux and Weidenmuller (1969), Philpott (1977)), which 
in the restricted space of configurations generated using the finite-depth potential, has 
been studied for the giant resonances and the radiative capture reactions probing the 
microscopic structure of these resonances (Barz et al (1977,1978), Fladt et al (1988)). 
This is insufficient for nuclei close to drip lines, where it is essential to have a most 
realistic description of bound state subspace. For that reason, the corner-stone of 
our approach, called the Shell Model Embedded in the Continuum (SMEC), is the 
realistic SM itself which is used to generate the A-particle wavef unctions. This choice 
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implies that the coupling between SM states and the one-particle scattering continuum 
must be given by the residual interaction. In our case, we use the residual interaction 
in the form : 

V = -V (a + bP? 2 )5(r 1 -r 2 ) , (1) 

with a + b = 1 and a = 0.73 . 

The key element of both SMEC and CSM is the treatment of single-particle 
resonances, which on one side may have an important amplitude inside a nucleus and, 
on other side, they exhibit asymptotic behaviour of scattering wavefunctions (Bartz 
et al 1977). The part of resonance for r < r c , where r c is the cut-off radius, 
is included in Q subspace, whereas the remaining part is left in the P subspace. 
The wavefunctions of both subspaces are then renormalized in order to ensure the 
orthogonality of wavefunctions in both subspaces. 

In the SMEC calculations, we solve identical equations as in the CSM (Bartz et 
al (1977)) but due to specificity of exotic nuclei, ingredients of these calculations are 
modified. For the bound states we solve the SM problem : Hqq^i — , using the 
code ANTOINE (Caurier (1989)). H QQ = QHQ is identified with the realistic SM 
Hamiltonian and are the A-particle (quasi-) bound wavefunctions. The quasi-bound 
resonances in the continuum are included as well. 

For the continuum part, we solve the coupled channel equations : 

(£(+) - H PP )i< +) = , (2) 

where index c denotes different channels and Hpp = PHP . The sign ± characterizes 
the boundary conditions, i.e., whether we consider incoming ' — ' or outgoing '+' 
scattering waves. In our case, we have ingoing wave in the input channel and outgoing 
waves in all channels. The structure of (A — 1) - nucleus is given by the SM, whereas 
one nucleon occupies a scattering state. The channel states are defined by coupling one 
nucleon in the continuum to a 'hole state' of (A — 1) - nucleus. The SM wavefunction 
has an incorrect asymptotic behaviour for states in the scattering continuum. On 
the other hand, the standard approach consisting of adjusting values of the two-body 
matrix elements to the experimental data and/or modifying the monopole term of 
the effective interaction (Dufour and Zuker (1996)), makes the definition of average 
single-particle field and the radial dependence of the single-particle wavefunctions 
somewhat arbitrary. Therefore, to generate both single-particle resonances and the 
radial formfactors of occupied orbits entering the coupling matrix elements (|l|) between 
states in the subspaces Q and P , we use the finite-depth average potential of Saxon- 
Woods type with spin-orbit part included. The parameters of the average potential 
are fitted to reproduce experimental single-particle states, whenever their identification 
is possible. This is the initial choice of potential, because the microscopic coupling 
of bound and scattering states generates the correction term to the mean-field which 
modifies the single-particle wavefunctions for each J n of the quasi-bound state . This 
correction term must be carefully taken into account otherwise, the solutions including 
coupling of Q and P subspaces could be non-orthogonal. 
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Figure 1. SM with Cohen-Kurath interaction and SMEC (labeled by residual 
interaction strength Vb) vs. experimental T = 1 states of 8 B nucleus. The proton 
threshold energy is adjusted to reproduce position of the ground state. The shaded 
regions represent the width of resonance states. 

The third system of equations are the coupled channel equations : 

(E+ - H pp )ujI +) = HpQ^i = w t (3) 

with the source term Wi , which is given by the SM structure of A - particle wavefunction 
for state <E>j . These equations define functions oj\ + ^ which describe the decay of quasi- 
bound state $j in the continuum. The source Wi couples the wavefunction of A - 
nucleon localized states with (A — 1) - nucleon localized states + one nucleon in the 
continuum. The formfactor of the source term is given by the same average potential 
as used in the P-space (||). The residual coupling ([!]) is also identical. The complete 
solution can be expressed by means of three functions : $j , and oji : 

n =&+ + - H^r 1 < i h qp i c E > , (4) 

where Hqq = Hqq + HgpGp H P q is the effective SM Hamiltonian including the 
coupling to the continuum and Gp is the Green function for the motion of single 
particle in P . Matrix Hqq is symmetric and complex. It can be diagonalized by the 
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orthogonal transformation : — > $j = J^fiji^i with complex eigenvalues E^ — \iTi. 
The eigenvalues of Hqq at energies E = Ei determine the energies and widths of 
resonance states. With this changes, one obtains : 

* c e = Ce + E^[£ ~Ei + (ifflfi]- 1 < $i | H | & > , (5) 

i 

for the wavefunction of continuum state modified by the discrete states and : 

fii = ^ + E/ dE£ B Ef (EW-E)- x <£ B Ef \H\Q i > , (6) 
for the wavefunction of discrete state modified by the continuum states. 



3. Results 

3.1. Spectroscopy of 8 B nucleus 

As an example of SMEC calculations for bound and resonant states, let us consider the 
proton-rich nucleus 8 B . Fig. 1 compares SM spectrum for T = 1 states of 8 B calculated 
in the p - shell using the Cohen and Kurath (1965) (CK) interaction, with those obtained 
in the SMEC for two different strength Vq of the residual interaction ([J). In the SMEC 
calculations, the Saxon- Woods average potential with parameters : Uq = —31.924 MeV, 
U so = 6.86 MeV, Rq = R so = 2.95 fm and a = a so = 0.52 was used as a starting guess. 
The spectrum of 8 B does not show large sensitivity to the choice of these parameters. 
This initial average potential is then modified for each J n state by the coupling to the 
continuum. This correction make the average potential deeper, it increases slightly the 
diffusivity and produceses a maximum at the center of the potential. It is the corrected 
potential which is then used to calculate the radial formfactors of coupling matrix 
elements and single-particle wavefunctions. The coupling to the continuum mixes the 
unperturbed SM states. For example for Vq = 1200 MeV-fm 3 in ([]]) , the wavefunction 
for the first if resonance has an overlap of 0.992, 0.1 + 0.005 i and 0.07 + 0.003 i with 
unperturbed SM states if, if and I3 , respectively. The calculated width of if state 
is 17 keV and 29 keV for Vq = 650 and 1200 MeV-fm 3 respectively. The experimental 
width for this state is 37 ± 5 keV. One should notice that also the energy difference 
between 2f and if states is improved by inclusion of the coupling to the continuum. 



3.2. Radiative capture reaction 7 Be(p, r y) 8 B at low energies 

The /3 + decay of 8 B, which is formed by the reaction 7 Be(p, 7) 8 B at the center of mass 
(cm.) energy of about 20 keV, is the main source of high energy solar neutrinos. In the 
absence of data in this region, the input of standard model for the solar neutrino problem 
(Bahcall and Ulrich (1988)) should be compared with those of various calculations. In 
the SMEC, the initial wavefunction \l/j([ 7 Be + p\j^ =1 ^ 
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Figure 2. Multipole contributions to the total capture cross section (left hand side) 
and the astrophysical S-factor (right hand side) of 7 Bc(p, 7) 8 B as a function of the 
center of mass energy. The SMEC calculations has been done with the residual 
interaction strength V = 650 McV-fm 3 . 



and the final wavefunction \I//([ 8 B] jj =2 +) is: 
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I t and s denote the spin of target nucleus and incoming proton, respectively. Aj^^ f is 
the coefficient of fractional parentage and is the s.p. wavefunction in the many- 

particle state Jf . With the wavefunctions ^j(r) and */(r) , we calculate the transition 
amplitudes: 

T EC = C(E£)i l "j f l b j b j a < C5J f m f | J imi >< l b 0£0 | l a > 
for El and E2 and : 

rpMl _ ila /jLN j f < lSJfmf I JiUli > 
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for Ml transitions, respectively. In the above formula, 5 = nii — nif, a = \/2a + 1, a = 
^a(a + 1) and l(^f*j b j b = J u i b j b rC 4 ) i^j a dr ■ The radiative capture cross section can then 
be expressed in terms of those amplitudes as : 
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Figure 3. Multipole contributions to the total capture cross section (left hand side) 
and the astrophysical 5-factor (right hand side) of 7 Bc(p, 7) 8 B as a function of the 
center of mass energy. The SMEC calculations has been done with the residual 
interaction strength V = 1200 MeV-fm 3 . 



[i stands for the reduced mass of the system. Fig. 2 shows the calculated multipole 
contributions to the total capture cross section and the astrophysical S-factor as a 
function of the cm. energy. The calculation is done for Vq = 650 MeV-fm 3 . The proton 
threshold energy is adjusted to agree energies of calculated and experimental 1 + state. 
The photon energy is given by the difference of cm. energy of [ 7 Be + p] j. =1 + system and 
the experimental energy of the 2f ground state of 8 B . The value of the S-factor at cm. 
energy of 20 keV: S = 22.58 eV-b , agrees with certain older experiments (Parker (1966), 
Kavanagh et al (1969)) but disagrees with the most recent one of Hammache et al (1998). 
The ratios of SMEC S factors at different cm. energies are: S(20)/S(100) = 1.07 and 
S(20)/S(500) = 1.05. The ratio of Ml and El contribution is: a M1 /a E1 = 8.13-10- 6 , 
4.55-10 -5 and 4.3-10~ 3 at 20, 100 and 500 keV, respectively. The resonant part of 
Ml transitions, which yields the contribution of S Ml = 24.72 eV-b at the if resonance 
energy, decreases fast and becomes S M1 = 9.6-10 -2 , 9.6-10 -4 , 1.8-10 -4 eV-b at cm. 
energies 500, 100, 20 keV, respectively. This value is proportional to the square of 
spectroscopic amplitude of p-states, which for the CK interaction is —0.352 and 0.567 
for pi/2 and p 3 / 2 respectively. Similar small values of spectroscopic amplitudes are 
obtained for Kumar (1974) and PTBME et al (1992) interactions (see Brown et al 
(1996)). At the position of 1 + resonance, the predicted by SMEC calculations value of 
S'-factor (S = 47.05 eV-b) is smaller than seen in some experiments (Kavanagh et al 
(1969), Vaughn et al (1970), Filippone et al (1983)). Part of this discrepancy could be 
related to the absence of quenching factor in our calculations. 

Fig. 3 shows the same as Fig. 2 but for V = 1200 MeV-fm 3 . Now, the calculated 
value of the S factor at cm. energy of 20 keV is S = 22.70 eV-b, rather close to 
the value obtained for Vq = 650 MeV-fm 3 , and disagrees with the recent experimental 
value given by Hammache et al (1998) . One should stress however, that different 
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experiments show significant fluctuations in the extrapolated values of S . The ratios 
of S factors at different cm. energies: S(20)/S(100) = f .07 and S(20)/S(500) = 0.997, 
show much stronger dependence on energy than seen for Vo = 650 MeV-fm 3 . The 
resonant part of Ml transitions at the if resonance energy, yields a somewhat smaller 
value : S M1 = 18.23 eV-b, than for Vo = 650 MeV-fm 3 . Finally, at the position of 
1 + resonance, the predicted by SMEC calculations value of S-factor (S = 41.97eV-b) is 
somewhat smaller than seen in the experiments and for V = 650 MeV-fm 3 . 

4. Conclusions 

In this work we have shown results of first calculations using the SMEC which couples 
the realistic SM solutions for (quasi-) bound states with the scattering solutions of 
one-particle continuum. The application to 7 Be(p, 7) 8 B reaction yields satisfactory 
description of different components of the radiative capture cross section, including 
the resonant components. In future, more unstable nuclei should be studied in the 
SMEC approach to systematically address the problem of effective interactions in the 
extreme conditions of exotic nuclei. At present, we are applying the SMEC to many 
other reactions of astrophysical interest such as 14 C(n, 7) 15 C, 16 0(p, 7) 17 F, 18 0(n, 7) 19 0. 

Acknowledgments 

Authors would especially like to thank S. Drozdz and I. Rotter for many clarifying 
discussions. We are grateful to E. Caurier and A. Lefebvre for encouragement and 
useful discussions. The work was partly supported by KBN Grant No. 2 P03 B 14010 
and the Grant No. 6044 of the French - Polish Cooperation. 



A study of nuclei of astrophysical interest in the continuum shell model 



9 



References 

Bahcall J N, Neutrino Astrophysics (Cambridge University Press, 1989). 
Bahcall J N and Ulrich R K 1988 Rev. Mod. Phys. 60 297 
Balashov V V et al 1964 Nucl. Phys. 59 414 
Bartz RW et al 1977 Nucl. Phys. A 275 111 

— 1977 Nucl. Phys. A 307 285 
Bcrtulani C A 1996 Z. Phys. A 356 293 
Brown B A et al 1996 Nucl. Phys. A 597 66 
Caurier E 1989 unpublished 

Cohen S and Kurath D (1965) Nucl. Phys. A 73 1 
Dobaczcwski J et al 1996 Phys. Rev. C 53 2809 
Dufour M and Zukcr A P 1996 Phys. Rev. C 54 1641 
Fano U 1961 Phys. Rev. 124 1866 
Filippone B W et al 1983 Phys. Rev. Lett. 50 452 
Fladt B et al 1988 Ann. Phys., NY 184 254, 300 
Hammachc F 1998 Phys. Rev. Lett, in print 
Julies R E et al 1992 S. Afr. J. Phys. 15 35 
Kavanagh W 1969 Bull. Am. Phys. Soc. 14 1209 
Kumar N 1974 Nucl. Phys. A 235 221 

Mahaux C and Weidenmuller H 1969 Shell-Model Approach to Nuclear Reactions (Amsterdam: North- 
Holland) 

Parker P D 1966 Phys. Rev. 150 851 
Philpott R J 1977 Nucl. Phys. A289 109 
Vaughn F J et al 1970 Phys. Rev. C 2 1657 



